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In this paper, we investigate the quest for a single target, that remains fixed in a lattice, by a set 
of independent walkers. The target exhibits a fluctuating behavior between trap and ordinary site of 
the lattice, whereas the walkers perform an intermittent kind of search strategy. Our searchers carry 
out their movements in one of two states between which they switch randomly. One of these states 
(the exploratory phase) is a symmetric nearest neighbor random walk and the other state (relocating 
phase) is a symmetric next-nearest neighbor random walk. By using the multistate continuous-time 
random-walk approach we are able to show that for dynamical targets, the intermittent strategy 
(despite the simplicity of the kinetics chosen for searching) improves detection, in comparison to 
displacements in a single state. We have obtained analytic results, that can be numerically evaluated, 
for the Survival Probability and for the Lifetime of the target. Thus, we have studied the dependence 
of these quantities both in terms of the transition probability that describes the dynamics of the 
target and in terms of the parameter that characterizes the walkers' intermittency. In addition to 
our analytical approach, we have implemented Monte Carlo simulations, finding excellent agreement 
between the theoretical-numerical results and simulations. 
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I. INTRODUCTION 

Search problems have recently experienced a rapid 
growth and motivated a great deal of work in the most 
various situations (see Ref. fl*! and references therein): 
fishermen and shoals, prey and predators, two molecules 
in the course of a reaction, a protein that looks for an 
specific site on DNA strand, medical drugs and illnesses. 
Thus, search problems spans over a wide range of do- 
mains and fields. When the target if fixed at a given 
location in space, the search problem is equivalent to the 
trapping problem, i.e., the situation where a set of walk- 
ers independently diffuse in space until one of them is 
caught by the trap. 

Among different forms of search strategies [11-01 , the so 
called intermittent strategies, which combine a phase of 
relocation (where the searcher may or may not be capable 
to capture the target), with a phase of search (where tar- 
get capture is always allowed), have been proved relevant 
and able to be optimized at various scales. Intermittent 
motion occurs in a wide array of living organisms from 
protozoans to mammals. It has been observed that nu- 
merous animal species switch between two distinct types 
of behavior while foraging or searching for shelter, or 
mate 0-01 • At a microscopic scale, we find intermittent 
motion, e.g., in the binding of a protein to specific sites on 
DNA for regulating transcription, as it is the case when 
the protein has the ability of diffuse in one dimension by 
sliding along the length of the DNA, in addition to their 
diffusion in bulk solution 

In Refs. [TT| - [l3j a theoretical model for the search ki- 
netics of a hidden target was presented, assuming that 
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each searcher could be in either of two states of diffu- 
sion. It was shown that intermittent strategies always 
improve target detection in comparison with the single- 
state displacement. An important aspect in the search- 
ing dynamics that has recently been studied [l^, is 
to consider not only the different states in the motion of 
walkers, but also what happens with the trapping process 
when additional dynamical effects are taken into account. 
For instance, in Ref. [l^, 13 the sighting range and the 
smell capacity of predators was considered as a sort of 
additional search ability. 

The aim of this study is to complete and to extend 
previous results [TT| - [l3j and to present another relevant 
case for modeling real search problems: The inclusion 
of a ffuctuating behavior in the target. These fluctua- 
tions can modify and prevent the encounter between a 
searcher and the target to be successful. This behavior 
may be due to the internal evolution of the trap or due 
to their interaction with a changing environment. For in- 
stance, in a chemical context, the activation or deactiva- 
tion of a reagent can be caused by external factors (pho- 
tons, solvent molecules, etc.) [l^. In biological contexts, 
the dynamical behavior of the target is also a determi- 
nant, e.g., reactions occurring within biomembranes re- 
quire some geometrical configurations in the biomolecule 
structure to be completed. The absent of these configu- 
rations inhibit the reaction, whereas stochastic changes 
in the molecule geometry can let it take place. Even the 
delivery of drug in medical treatments can involve block- 
ing chemical reactions, in order to boost the delivered 
medicine effectiveness [iq . 

Ref. [Trj introduced a generalization of the trapping 
model which allows encounters between particles of two 
kinds, A and B, with or without annihilation depending 
on the internal state of the particle A. Particle A, which 
is identified as the trap in this work, has two states: An 
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active one in which the annihilation does take place, and 
an inactive one in which it behaves as a regular site. The 
particles B are our walkers. In this work, we take a step 
forward in modeling search problems by the formulation 
of an unified framework which comprises the dynamical 
behavior of the trap and the intermittent search strat- 
egy performed by the walkers. We exploit the theory of 
multi-state random-walk (RW) ^IS] , we use the concepts 
of Survival Probability (SP) and Mean Target Lifetime 
(MTL), and we establish the connection to the First- 
Passage Time (FPT) corresponding to the problem of 
one walker. 

The outline of this paper is as follows. The next sec- 
tion presents our model and gives the basic definitions 
and concepts. Also, this Section describes the analytical 
approach to the trapping process. Section III presents 
the main results for the SP and the MTL of the target 
through a comparison between the numerical evaluation 
of our analytical framework and Monte Carlo simula- 
tions. In Sec. IV we discuss our conclusions. Finally, 
in Appendix A we develop the analytical calculations 
of Sec. Ill, corresponding to infinite chains and rings, 
whereas in Appendix B we consider the high transition 
regime for the trap. 



II. ANALYTICAL APPROACH 

A. The Model 

We restrict our work to chains (finite and infinite) and 
assume that the dynamical trap is held fixed at the origin 
of the lattice. A set of walkers, uniformly distributed 
along the chain, starts the "search" at i = 0. At the trap 
site, the following situations may occur: 

• The trap is in an active status, i.e., it works (the 
first walker reaching the trap is caught with prob- 
ability one, i.e., perfect trapping). 

• The trap is in a passive status and stays that way 
until the searcher leaves it, i.e., the trap behaves 
like any other chain site, and capture can not be 
carried out. 

• The trap is in a passive status but changes its con- 
dition before the searcher leaves, i.e., the capture 
is also performed. 

We denote the internal states of the trap by i = 1 (active 
status) and i — 2 (passive status). On the other hand, we 
assume that each predator can make two types of motion 
on the lattice: 

• Exploration: RW with symmetric jumps to first 
nearest neighbors, with transition probability per 
unit time A, and 

• Relocalization: RW with symmetric jumps to sec- 
ond nearest neighbors, also with transition proba- 
bility per unit time A. 



We also assume that the walkers' dynamics and the dy- 
namical behavior of the trap are independent. 

The proposed composite process can be described by 
the coupled master equations 

5Fi,»o(s,i|so,0) 



dt 



^P2.^„is,tWo,0) 

dt 



= APi,,„ (s, t|so, 0) + 72 P2A0 (s, t\sQ, 0) 
-7iPi^,„(s,t|so,0), (1) 

AP2,*o {s, 1 1 So, 0) + 71 Pi.io (s, t|so, 0) 

-l2P2,^o{s,t\so,0), (2) 

where Pi^ig{s,t\sQ,t — 0) is the the conditional probabil- 
ity of the walker of being at site s with the trap in state 
i at time t, given that it was at site sq with the trap in 
state io at < = 0. For simplicity, we have restricted the 
activation - deactivation process of the trap to time ex- 
ponential density functions with parameters 7^, i.e., ji is 
the probability transition rate of the trap to make a tran- 
sition from its state i to the other state. The dynamical 
evolution of the walkers, taking into account its intermit- 
tency, is described by the operator A. Particularly, for a 
chain, we get 

[A],,,' = ^ [(1 - a)iS,^s'-i + Ss,s'+i) 



(3) 



where a is the parameter that regulates the walker's fre- 
quency intermittency and A its diffusion constant (see 
Fig. p. 





FIG. 1. (a) Schematic transitions of the walker to/from site 
s (away from the trap: s 7^ —1,0,1) and (b) Walker tran- 
sitions to/from s = (trap site). A walker dwelling at site 
could be trapped with rate 72 (the probability transition 
rate for activation of the trap). The dynamics of the trap is 
independent of the dynamics of walkers. 



B. The Trapping Process 

We will focus on the SP of the dynamical target, i.e., 
the probability that the target remains undetected up to 
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a time t, and its closely related quantity, the MTL [i9| . 
which compute the time in which the first walker reaches 
the target under the appropriate circumstances of cap- 
ture. We define Fi^^^ (0, i[so, 0) as the First Passage time 
density through the site at time t, when capture is 
possible, given that the searcher was at sq with the tar- 
get in state iq at time t ~ 0. In the way of Ref. Ig*], 
we introduce the notion of generalized state which takes 
into account the position of the walker and the state of 
the target, {s,i). The connection between FPT density 
at (0,1) at time t from (so,jo), Fi i„{0,t\sQ,0), and the 
conditional probability Png{s,t\sQ,t — 0) is established 
by 



The MTL is defined in finite domains as [l^ 



Tn= / <i>N{t)dt, 
Jo 



and in the bulk limit as 



Tp- / $p(t)dt. 
Jo 



(8) 



(9) 



We left for the Appendixes the detailed calculations of 
the magnitudes presented in this section for the cases of 
infinite chains and rings of M sites. 



A,i(0,w|0,0) 



(4) 



which is the known Siegert's formula [20|, generalized to 
internal states '21']. We are denoting the Laplace trans- 
form of a function of i by a caret over the corresponding 
function. Thus, for example, 

p2,iois,u\so,0) = C{Pt,i„{s,t\so,0)} 

e-"*P,,,„(s,t|so,0)dt. 



When trapping occurs independently of the initial 
state of the target, the SP in presence of only one walker 
may be written (if so ^ 0) as 



$i(0,t|so,0) = l-^( 

40 = 1 



i^l,.„(0,T|so,0)dT. (5) 



6ig is the probability of the initial state of the target. 
Thus, the target is initially active with probability 61 or 
inactive with probability 02- The SP at time t, $Ar(i), 
of the dynamical target at the origin in the presence of 
N independent walkers that diffuse on an M-sites lattice 
can be written in terms of th(3 SP in the presence of only 
one walker, $i(0, i|so,0) as [19|] 



N 



$iv(0 = 1 - 



Af - 1 



^(l-<i>i(0,t|so,0)) 



(6) 



where we have assumed a uniform probability distribu- 
tion for the initial position of the walkers, i.e, the prob- 
ability that a given walker is initially at a particular site 
50(7^ 0) is {M — Notice that we explicitly exclude 

the possibility of having a walker at the position of the 
target at t = 0. In the bulk limit, N ^ 00, M ^ 00, with 
N/M — !> p, where the constant p is the concentration of 
walkers, we get 



$p(i) =exp -p^(l-$i(0,t|so,0)) 



(7) 



III. RESULTS 

In this section, we illustrate the general framework 
presented in the previous section. We consider one- 
dimensional systems and give some general ideas to in- 
terpret the obtained results. The inverse Laplace trans- 
form involved in the analytical expressions, given in the 
Appendixes, is evaluated numerically [1^ for obtaining 
concrete results and then we establish a comparison with 
independent Monte Carlo simulations. 

A brief review of our simulation methodology is 
appropriate at this point. We uniformly distribute 
the searchers (with probability p per site) in a one- 
dimensional lattice with periodic boundary conditions. 
The target is placed at the origin of the lattice. The 
propagation of the searchers in the presence of a dy- 
namical target is implemented as follows. Each searcher 
has assigned an internal clock (all start synchronized at 
time t = 0) which is updated according to their waiting 
time probability distributions. For the activation - deac- 
tivation process of the target a similar procedure to the 
searchers is used; the target has assigned his own inter- 
nal clock which is updated with time exponential density 
functions with parameters 71 and 72. We define an in- 
dicator function that records the needed information: if 
the target was captured up to a certain time (for the 
SP) and whether the target was captured and the time 
in which this happened (for the mean target lifetime). A 
randomly chosen walker take a step, to its nearest neigh- 
bors with probability (1 — a) or next- nearest neighbors 
with probability a and left or right with equal probability 
(1/2). We check if the trapping conditions are fulfilled, 
and if it does, we stop the dynamics, update our indica- 
tor function, and generate a new ensemble of walkers. If 
it does not, we continue the dynamics by taking another 
randomly chosen walker. Again, if trapping occurs, the 
indicator function is updated and the dynamics stopped; 
if not, the walk continues. The output of interest of each 
realization is, for the SP, whether it was captured up to 
a certain (predefined) time, and the time of capture for 
the mean target lifetime. 

In this section, we show the numerical results obtained 
from the analytical expressions for the infinite chain (see 
Appendix lA 1[) and the finite ring (see Appendix lA 21) . 
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FIG. 2. Analytical-Numerical calculations (lines) and Monte 
Carlo simulations (symbols) for the SP, <l?jv(a;t), up to time 
t — 20 for different target transition rates (7). (I) (diamonds) 
7 = 0.01, (II) (triangles) 7 = 0.1, (III) (circles) 7 = 1 and 
(IV) (squares) 7 = 10. We have also included for comparison 
the static trap case (thick solid line). 
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FIG. 3. Analytical-Numerical calculations (lines) and Monte 
Carlo simulations (symbols) for the Mean target lifetime 
(MTL) Tp, for different target transition rates (7). (I) (di- 
amonds) 7 — 0.01, (II) (triangles) 7 — 0.1, (III) (circles) 
7 = 1 and (IV) (squares) 7 — 10. We have also included for 
comparison the static trap case (thick solid line). 



For both the infinite and finite cases a searchers' concen- 
tration p ^ 0.1 was used (for the finite case the concen- 
tration is defined by p = N/M). All the finite chains 
considered in the figures correspond to a ring of M = 20 
sites, with exception of Fig. |4l where different sizes (M) 
of the ring are explicitly stated. All times are given in 
units of the inverse of the diffusion constant (1/A). It is 
worth to comment that when we talk about target transi- 
tion rates 7^, these are low (high) relative to the diffusion 
constant A. An equivalent interpretation can be made if 
we consider the target mean sojourn time in state i as 
7~^; this will be long (short) on the time scale determined 
by the propagator A (A^^). In the following we will 
consider symmetric transition rates for the activation- 
deactivation process of the target, 71 = 72 = 7- 

In Fig. [2] we present curves (for the finite case) corre- 
sponding to $jv(a,t) for a fixed evolution time t — 20. 
Notice how the intermittent search can improve the de- 
tection probability, i.e, minimize the SP of the target, 
compared with the single state search {a 0, a 1). 
As a comparison we also include the "static trap" case, 
i.e., the target is always active. As can be seen from the 
figure, an optimal value for a can be found for each target 
transition rates 7 chosen. Even though all curves present 
a similar behavior, it is apparent that the transition rate 
7 plays an important role. The ratio between the max- 
imum value of the SP (at a = 1) and it's minimum is 
almost of 80% (120%) for 7 = 0.1 (7 = 0.01). At high 
values of 7, the "static trap" case is approached. 



The curves shown in Fig. [3] correspond to the MTL, in 
the finite function of the walker intermittency 

parameter a, for different target transition rates 7. As it 
is clear from the figure, the results show the same trend 
as the SP (Fig. [2]) revealing also a remarkable rise in MTL 
for low values of the target transition rates 7. It could 
be inferred from the curves that with a modest transition 
rate value (7 — 0.1) the target could almost double its 
lifetime expectancy while a high "activity" of the target 
(7 > 1) leads it to the static case. Note that although 
MTL has less information than the SP, it shows to be a 
simple and efficient tool for characterizing the proposed 
search scheme. 

Figure [3] depicts the behavior of the MTL for a fixed 
target transition rate (7 = 1) as a function of the walker 
intermittency parameter a and for different sizes (M = 
20, 40, 60, 100, 200, 1000, 00) of the chain. In all cases was 
used a concentration of searchers p = 0.1. Notice how the 
finite chain (ring) approaches the infinite chain even for 
values of M not too large. As can be seen from the figure, 
the minimum in MTL is maintained for all system sizes, 
which constitutes a robust property of the intermittent 
search approach. In Fig. [5] we consider the behavior of 
the SP, ^]\[{a;t), in the high transition regime of the 
dynamic target for a fixed evolution time t — 20. In this 
limit, the behavior of the SP approaches an imperfect 
trap (see appendix [B|) . with i/ = (71 -I- 72)72/71 being a 
"measure of the imperfection" of the trapping process. 
When ^ there is no trapping and if — 00 perfect 
trapping is achieved. Notice how the dynamical trapping 
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FIG. 4. Analytical-Numerical calculations for the Mean Tar- 
get Lifetime (MTL) for a fixed target transition rate (7 = 1) 
as a function of the walker intermittency parameter a and 
for different sizes M of the chain. From bottom to top 
M = 20, 40, 60, 100, 200, 1000 and M = cx) (thick solid hne). 



resembles the imperfect case even for values of 7^ not 
too large. It is worth remarking the excellent agreement 
between the analytical-numerical results and the Monte 
Carlo simulations. 
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FIG. 5. Analytical - Numerical calculations (lines) and Monte 
Carlo simulations (symbols) for the SP, '^(a.\t\ up to time 
t — 20, in the high transition regime of the target. Lines with 
symbols corresponds to the 'imperfect' (high transition) case 
and the same type of line (without symbols) for the dynamical 
case. For instance (I) (squares) f = 2, 7 = 1; (II) (circles) 
= 5, 7 = 2.5 and (III) (squares) 1/ — 20, 7 = 10. We have 
also included for comparison the static case (thick solid line) 



IV. CONCLUSIONS 

We have presented a simple model for the search ki- 
netics of a set of walkers performing intermittent motion 
in quest for a dynamical target. The model is based only 
on RW, and our results complement and extend previous 
related results given in Ref. [lll - [l3| . However, this model 
differs from those mentioned. In our previous work, the 
first searcher that finds the target, captures it with proba- 
bility one (we denominate that situation the "static trap" 
in this work). In the present work, an encounter walker- 
target does not necessarily end in capture, but depends 
rather on the state of the dynamic target. 

We have considered the target's survival probability 
(at a fixed time) and the target's lifetime, and also stud- 
ied the dependence of these quantities on both the tran- 
sition probability (7^) between the states of the target 
and the parameter that characterizes the walker's inter- 
mittency (a). Thus, we have established that the SP is a 
non-monotonic function of a for a wide range of the tran- 
sitions probabilities 71 and 72, showing that intermittent 
strategies still improve target detection when compared 
with the single-state displacement. This confirms the 
utility of the intermittent search approach [2^ even in 
the case of a dynamical target. 



We introduced the MTL and its connection with the 
SP was established. As it was the case for SP, MTL was 
also a non-monotonic function of a for several values of 
the transitions probabilities 7^, adequately depicting the 
improvement provided by the intermittent search strat- 
egy. Although MTL carries less information than the 
SP, it has shown to be an efficient global optimizer for 
search strategies using intermittent motion. In all cases 
the agreement between analytical-numerical results and 
Monte Carlo simulations was quite good. 

Thus, we have fulfilled our goal of presenting a simple 
model, based only in diffusion, that captures in an unified 
framework the dynamical behavior of the target and the 
intermittent search strategy performed by the walkers. 
The present scheme is both simple enough to be studied 
analytically and rich enough to be able to mimic the in- 
fluence of the target's dynamics in the capture process 
and it shows that intermittency is always favorable for 
optimizing the search. 

The present model of intermittent search can be gen- 
eralized in several directions: higher dimensions, contin- 
uous systems, non-Markovian target dynamics, etc. All 
of these aspects will be the subject of future work. 
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Appendix A 

Here, we describe some essential details of the calcula- 
tions of Section (jll B|) . We focus on the Pi^i„{s,t\so,t = 
0), which are the building blocks for the SP and MTL. 
Given that we now apply our model only to chains, we 
drop out the vector notation. 



1. Infinite Chain 

The solution of Eq. ([Ij, with the particularization of 
Eq. for an infinite (homogeneous) chain, can be 
given following the guidelines of Refs. [13, |2i|. The 
exploitation of the indicated formalism leads us to an- 
alytic closed expressions for Piig{s,t\sQ,t = 0) in the 
Fourier-Laplace space. However, as we are interested in 
Pi^ig{s,t\so,t — 0) (we take i = 1 as the active status of 
the target) we only show the solution for this case. The 
relevant transformed results reads 

Pi,i(fc, u) = i {l2p°ik, u) + 7i^°(fc, u + T)), 

PiMk, u)^^ {l2P\k, u) ^ 72P°(fc, U + T)), (Al) 

where F = 71 + 72, and P^{k,u) = l/(u — A{k)), is 
the Fourier-Laplace transform of the conditional prob- 
ability P°{s,t\so,t = 0), corresponding to the inter- 
mittent walker being at site s at time i, given that it 
was at site sq at t = (without trap) and A{k) — 
A[(l — a)cosfc + acos2fc — 1] is the Fourier transform 
of the evolution operator A given by Eq. ([3]) . 

We are interested in obtaining results for any initial 
state of the target/trap. Therefore is useful to evaluate 
the average 

A,^o (fc, U) = Pis{k, U) 01 + A,2(fc, U) 02 

io 

- i (72P"(fc, u) + (7101 - ^202)P°{k, u + F)) ,(A2) 

where 0i is the target probability of being initially active, 
and O2 of being inactive and satisfy 0i + 6*2 = 1. As 
usual, we choose for 6i the equilibrium probabilities [31 ; 
di = 72(71 + l2)~\ O2 = 71(71 + 72)""^- The Fourier 
inversion of Eq. (|A2p could be calculated in an exact way 
resulting 



io 



1,^0 



(s,m|so, 0)6*^0 = G 



\/x\ - 1 ^ 



where 771 = xi - ^Jx\-\, 772 = 2^2 + \Jx\-^, G 
72/2AFQ!(a;i — ^2) and 



1-a , 1 \(\-a\ u + \[\ + a) 
^^'^ = -^^2Vl^j ■ 

Averaging Eq. (IA3|) over the starting positions (uniformly 
distributed) of the walker we arrive at 



E 

io = l,2 

50 5^0 



Pl,io(s,M|so, 0)61,^0 

G 2771 G 2772 



(A4) 



- 1 1 - m ^Jx\ - 1 1 - r\2 

Taking s = 0, so = 0, and considering Eqs. (|A3p 
and (IA4I). we can write 



^/:{(i-<i>i(o,i|so,o))}-- A,.o(o,w|so,o)a, 



1 



1 



40 = 1,2 
so#0 



^^Pl,l(0,w|0,0) ^^tf,2 

(A5) 

Eq. (jASp constitutes one of our main results and it allows 
us derive SP from Eq. ^ and MTL from Eq. ©, after 
taking the inverse Laplace transform. However, despite 
being able to obtain analytical results in Laplace space 
for Eq. (jA5l) , its length and complexity made the analysis 
a difficult task. The analytical inversion of the Laplace 
transform of the results seems to be beyond our possibil- 
ities, so we have used a numerical procedure [22] for its 
calculation in Sec. IHIl 



2. Ring of M sites 

For the finite case, we take the results from the pre- 
vious section, and obtain for the ring a solution in the 
form fH 



A^4(s,M|sO,i = 0) - Y P^,^ois + lM,u\so,t = 0). 

/ — — oo 

(A6) 

In order to evaluate the sum proposed in Eq. (jA6p . we 
focus our attention in one term of Eq. (|A3p . Given that 
for alM 0, \l\ M > (s-sq) and if / < \{s-so) + lM\ = 
— (s ~ So) — I M, we get 



Y 77l^-^°+"'l = 77l^-^"l + (77il^-^''l +77^'^"^"') Y 

Js-so| ^ ^M-\s-so\ ^ (A7) 



1 = 1 



l-77f 



M 



(A3) 



Working in a similar way with the other terms, we obtain 
the complete solution of Eq. (jA6p , for the state of capture 
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{i = 1) as 

oo 



G 



' |s-so| , M-|s-so 
'/I "I" '/I 



s — Sol I M— Is — Sol ' 



(A8) 



J 



For a uniform distribution of walkers along the ring is useful the expression 

2 



^ «ko, 0)0.0 =G 



?72 - 111' 



O(i-m) v^f^(i-r?r)(i-r?2) 



So5^0 



(A9) 



With Eq. (jA8l) valuated in s = 0, sq = 0, and taking into 
account (IA9I) we can write 



^£{(1- $1(0,^150,0))} = i ^ Fi'!!o(O-"ko,O)0,„ 



so/O 



jo = l,2 
soT^O 



1 



1 



z^Pi^l(0,u|0,0),^^4r.2 

so#6 

(AlO) 



From Eq. (|JT0)) . the SP (Eq. ©), $Ar(t), and the MTL 
(Eq. dS])), Tat, are obtained. However, as in the previous 
section, the size and complexity of Eq. (|A10p makes the 
inversion of the Laplace transform beyond our possibil- 
ities, so we need to use a numerical procedure [22] for 
obtaining the concrete results presented in Sec. IIIIl 



Appendix B 

In this appendix we consider the high transition regime 
in dynamical trapping, i.e., the behavior of the SP in the 
limit r 3> A. The calculation may be carried out starting 
with the Laplace transform of Eq. ([SJ , 



3>i(0,u|so 



,O) = i|^l-^i^i,.o(O,7.|so,O)0.oj • (Bl) 



r 



Let us consider the second term on the right hand side 
of Eq. (HI]) 



A,.o(0, u\so, 0)9,, = 72 P\0, w|so, 0) / 

io 

(72P"(0,u|0,G) + 7iP°(0,M + r|0,0)) . (B2) 

In the considered limit (F > A) , P°(0,u + F|0,0) ~ 
1/F [13, then 



^A,.o(O,^i|so,O)0, 



72P°(0,^ko,0) 
72PO(0,M|0,0)+7i/r 

vP°{OMso,0) 
l + iyP"(0,u|0,0) ' 



(B3) 



where v — T^2 / 7i • Notice that Eq. (|B3p adequately pro- 
vides the limits of perfect trapping (72/71 ^ 1, i.e., v — >■ 
00) E.oA,.o(0,u|so, 0)0,0 = P0(0,w|so,0)/PO(0,M|0,0) 
and no target/trap present (72/71 -C 1, i.e., v 0) 
J2,„ K^o (0, u\so, 0) = 0. Using Eq. ^ and Eq. (EH) 
we finally obtain 



$i(0,u|so,0) 



i/pO(0,w|so,0) 
1-f j/PO(0,u|0,0) 



(B4) 



This result resembles the general case when detection 
of the target upon encounter is less than certain, i.e., 
imperfect trapping 3]. From Eq. (IB4[) and using the 
same procedure from Appendix (|A 2|) the SP (Eq. [6|), 
$Ar(i), and the MTL (Eq.E]), Pat, could be evaluated for 
the present regime. 
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